A comparison of cosmological models using recent supernova data 
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We study the expansion history of the universe up to a redshift of z=1.75 using the 194 recently 
published Snia data by Tonry et. al. and Barris et. al. In particular we find the best fit forms of 
several cosmological models and H{z) ansatze, determine the best fit values of their parameters and 
rank them according to increasing value of Xmin (the minimum value of for each H{z) ansatz). We 
use a prior of Qom = 0.3 and assume flat geometry of the universe. No prior assumptions are made 
about validity of energy conditions. The fitted models are fourteen and include SCDM, LCDM, 
dark energy with constant equation of state parameter w (quiessence) , third order polynomial for 
H{l + z), Chaplygin gas, Cardassian model, w{z) = wq + wiz, w{z) = wo + zwi/{l + z), an oscillating 
ansatz for H{z) etc. All these models with the exception of SCDM are consistent with the present 
data. However, the quality of the fit differs significantly among them and so do the predicted forms 
of 'w{z) and H{z) at best fit. The worst fit among the data-consistent models considered corresponds 
to the simplest model LCDM (Xmin = 198.7 for Qom = 0.34) while the best fit is achieved by the 
three parameter oscillating ansatz (Xmin = 194.1). Most of the best fit ansatze have an equation of 
state parameter w{z) that varies between w{z) ~ —1 for z < 0.5 to w{z) > for z > 1. This implies 
that the sign of the pressure of the dark energy may be alternating as the redshift increases. The 
goodness of fit of the oscillating H{z) ansatz lends further support to this possibility. 
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I. INTRODUCTION 



One of the fundamental goals of cosmology is the un- 
derstanding of the global history of the universe. Us- 
ing objects of approximately known absolute luminocity 
(standard candles) in the nearby universe provides the 
current rate of expansion. Using more distant standard 
candles like type la supernovae (SnIa) makes it possi- 
ble to start seeing the varied effects of the universe's 
expansion history. Such cosmological observations have 
indicatedPl that the universe undergoes accelerated ex- 
pansion during recent redshift times. This accelerating 
expansion has been attributed to a dark energy compo- 
nent with negative pressure which can induce repulsive 
gravity and thus cause accelerated expansion. The sim- 
plest and most obvious candidate for this dark energyf^ 
is the cosmological constant |3J with equation of state 

w = ^^-l. 

p 

The extremely fine tuned value of the cosmological con- 
stant required to induce the observed accelerated expan- 
sion has led to a variety of alternative models where the 
dark energy component varies with time. Many of these 
models make use of a homogeneous, time dependent min- 
imally coupled scalar field (/) (quintessence^, @) whose 
dynamics is determined by a specially designed potential 
V{(j)) inducing the appropriate time dependence of the 

field equation of state w{z) — Given the observed 

w{z), the quintessence potential can in principle be de- 
termined. Other physically motivated models predictin 
late accelerated expansion include modified gravityli" 



Ig, Chaplygin gas 9], Cardassian cosmol ogy pOj. theories 
with compactified extra dimensions |lll ll^. braneworld 
models 13] etc. Such cosmological models predict specific 
forms of the Hubble parameter H(z) as a function of red- 
shift z in terms of arbitrary parameters. These parame- 
ters are determined by fitting to the observed luminocity 
distance ^^(z) using the relationO [H [3 (vahd m a 
fiat universe) 



^W = c[A(Mf); 

dz 1 + z 



(1.1) 



This is easily derived using the relation between di^{z) 
and the comoving distance r{z) (where z is the redshift 
of light emission) 

dL{z) ^ r{z){l + z) (1.2) 

and the light ray geodesic equation in a flat universe 
c dt = a{z) dr{z) where a{z) is the scale factor. 

Another similar approach towards determining the ex- 
pansion history H[z) is to assume an arbitrary ansatz 
for H{z) which is not necessarily physically motivated 
(it is 'model independent') but is specially designed to 
give a good fit to the data for dL{z). Given a particular 
cosmological model (or ansatz) for H{z; oi, a„) where 
ai, ..-Tan are model parameters, the maximum likelihood 
technique can be used to determine the best fit values of 
parameters (with la — 2a errors) as well as the goodness 
of the fit of the ansatz to the data. This technique can be 
summarized as follows: The observational data consist of 
N apparent magnitudes mi{zi) and redshifts zt with their 
corresponding errors Srui and Szi. Each apparent magni- 
tude is related to the corresponding luminocity distance 
dL of the SnIa by 
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m{z) — M + 5 log 



10 [ 



,dLiz). 
Mpc 
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(1.3) 
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where M is the absolute magnitude which is assumed 
to be constant for standard candles like Type la Snia. 
From equations Hl.l|l and (jl.^ll it becomes clear that the 
luminocity distance d^^z) is the 'meeting point' between 
the observed apparent magnitude m{z) and the theoret- 
ical prediction H{z). 

The apparent magnitude can also be expressed in 
terms of the dimensionless 'Hubble-constant free' lu- 
minocity distance defined by 



Dl{z) 



HodLiz) 



(1.4) 



as 



m(z) = M{M, Ho) + 5 log,oiDL{z)) (1.5) 

where M is the magnitude zero point offset and depends 
on M and Hq as 



M^M + 5log,o{^^) + 25 
lAlpc 



(1.6) 



The zero point offset is an additional model independent 
parameter that needs to be fit along with the model pa- 
rameters oi, ...,a„. However, since M is model indepen- 
dent its value from a specific good fit can be used directly 
to other fits of model parameters. Thus the observed 
mi{zi) can be translated to D'^f{zi) using equation (|1.5() 
for the best fit value of Mobs obtained from nearby SnIa. 
The theoretically predicted value D*^{z) in the context 
of a given model H{z; ai, a„) can be obtained by in- 
tegrating the equation Hl.l(l as 



Dl^{z) = (l + z) / dz' 



Ho 



H{z';ai, ...an) 



(1.7) 



The best fit values for the parameters ai, a„ are found 
by minimizing the quantity 



N 



{logioDt'{zi)-logioDf{zi))^ 



(1-8) 

where cr^ is the la redshift uncertainty of the data 
and (Tiog-^gDi^izi) is the corresponding la error of 
logioDl'^iz,). \ 

The probability distribution for the parameters 
oi, a„ is 5^ 

P(ai,...,a„)=AAe-x'(«i-'^")/2 (1.9) 

where TV is a normalization constant. If prior informa- 
tion is known on some of the parameters ai, ...,a„ then 
we can either fix the known parameters using the prior 
information or 'marginalize' i.e. average the probability 
distribution (|1.9|) around the known value of the param- 
eters with an appropriate 'prior' probability distribution. 
Here we use the former approach (fix the parameters with 
prior information) for simplicity. This simplification has 
negligible effect on our results as it can be verified by 



comparing some of our results with corresponding re- 
sults in the literature where marginalization has been 
used (e.g. Ref. for LCDM). 

It is straightforward to minimize Y^ (qi , an) us- 
ing numerical libraries like NAG jl^ l20l (see also 
[131) or packages like Mathematica [21| to find 
Xmm(ai> •■•:an)|2a| where 

Xmin is the minimum obtained 
for the best fit parameter values ai, a„. If Xmin/ ~ 
n) ^ 1 the fit is good and the data are consistent with 
the considered model H{z; ai, a„). 

The variable x^j„ is random in the sense that it de- 
pends on the random data set used. Its probability dis- 
tribution is a distribution for N — n degrees of freedom 
[r^ . This implies that 68% of the random data sets will 
give a such that 

X^(ai, a„) - x^(ai, a™) < AxL(") (1-10) 

where l^~xi„{n) is 1 for n — 1, 2.3 for n = 2, 3.53 for 
n = 3 etc. Thus equation (|1.1U|) defines closed ellipsoidal 
surfaces around ai, a„ in the n dimensional parameter 
space. The corresponding la range of the parameter ai is 
the range of a^ for points contained within this ellipsoidal 
surface. Similarly, it can be shown that 95.4% of the 
random data sets will give a such that 



,)-x'(ai,-,an) < AxLH (l-H) 



where lS.x^^{n) is 4.0 for n — 1, 6.17 for n — 2, 8.02 for 
n = 3 etc. Thus equation (|l.ll|l defines the 2a ellipsoidal 
surface in parameter space and similarly for higher cr's. 



II. COSMIC EXPANSION HISTORY 

We now apply the above described maximum likeli- 
hood method using a recently published data set consist- 
ing of 194 [N = 194) SnIa ElEl- This is a subset of 
the total of 253 published SnIa sample obtained by im- 
posing constraints Ay < 0.5 (excluding high extinction) 
and z > 0.01 (reducing peculiar velocity effects). Each 
data point at redshift z includes the logarithm of the 
Hubble-free luminocity distance logio{cD'^'^ (z)) and the 
corresponding error crjogioDilz) ■ ^ table of the data we 
used can be downloaded in electronic form These 
Hubble-free luminocity distances are obtained assurning 
a best fit value for the zero point magnitude offset M [l^ . 
We adopt this same value for M and choose not to treat 
M as an additional free parameter to fit (and marginal- 
ize) along with the parameters of each theoretical model 
studied. In the Appendix we demonstrate that marginal- 
ization over M would have negligible effect (0(1%)) on 
our results. Also comparison of our results for LCDM 
and quiessence (w(z) = constant) with the correspond- 
ing results of Ref. [ig where marginalization of M was 
implemented indicates that our simplified approach has 
negligible effect on the obtained results. This same con- 
clusion has also been reached in Ref. [2^ and it's origin 
is demonstrated in the Appendix. 
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In the construction of using equation H1.8() we have 
used a value of corresponding to uncertainties due to 
pecuHar velocities with Aw = A(cz) — 500km/ sec which 
implies cr^ — Az = {500km/ sec) /c. The minimization of 
(|1.8() was implemented for each theoretical model using 
a simple Mathematica code which can be downloaded 
along with the table of the data set used 22] (or can be 
sent by e-mail upon request). 

We now proceed to apply likelihood testing to various 
theoretical models. Each model is defined by its pre- 
dicted Hubble-parameter H{z). For example for LCDM 
we have 

H\z;norn) = i^f = H^[nom{l + z)^ + {l-norn)] (2.1) 

and there is a single parameter f2om to be fit from the 
data. The simplest model to consider however is SCDM 
defined by 

H'iz)=H',il + z)' (2.2) 

with no free parameters. Using equation (|2.2|) in 1)1. 7|l we 
calculate D]!/{z). We may then find using equation 
(|1.8|) and minimize to find Xmin- the SCDM case there 
are no free parameters to vary and no minimization is 
needed. We thus find = Xmin = 431.4 which implies 
Xmin/dof = 2.2 (do/=degrees of freedom). Since this 
value X^in/dof is significantly larger than 1 we conclude 
that SCDM does not provide a good fit to the Snia data. 

The next simplest model consistent with the fiatness 
indicated by WMAPp^ is LCDM defined by equation 
(|2.1|) . It is straightforward to evaluate D*^{z) numeri- 
cally (using equation 1)2. l|l in (|1.7(l ^ and use it to evalu- 
ate x^(^Om) from equation (|1.8|) . A minimization of this 
expression leads to|2^ 

Xrmn = xH^om - 0.34) = 198.74 (2.3) 

which implies x'^/dof = 1.03 {dof = 194-1 = 193). This 
model is clearly consistent with the data since /dof ~ 
1. The ItT errors on the predicted value of fiom = 0.34 
are found by solving the equation 

xH^mla)~xLn^^xUn^l)^l (2.4) 

which leads to 

f2om = 0.34 ± 0.032 (2.5) 

This result is identical with the result of Ref.^3 even 
though our la errors are slightly smaller. We note here 
for comparison with the models discussed below (where 
the prior flom — 0.3 is used) that x'^i^Om = 0.3) = 200.3. 

In Fig. 1 we show a comparison of the observed 194 
SnIa Hubble free luminocity distances along with the the- 
oretically predicted curves in the context of SCDM (con- 
tinuous line) and LCDM (dashed line). In this case it is 
even visually obvious that LCDM provides a good fit to 
the data contrary to the case of SCDM. This visual dis- 
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FIG. 1: The observed 194 SnIa Hubble free luminocity dis- 
tances along with the theoretically predicted curves in the 
context of SCDM (continuous hne) and LCDM (dashed line). 

tinction is not possible when comparing the other data- 
compatible models discussed below with LCDM. We thus 
do not attempt to include the theoretical curves corre- 
sponding to other models on the same plot. 

We now consider other more general models and 
ansatze which however reduce for certain parameter 
values to LCDM. If these parameter values give a 
X^{LCDM) that is beyond the 2a level away from the 
minimum y^,„ then we would conclude that LCDM is 
disfavored compared to the better fit model. Even if we 
just find models with x^in < xliini^CDM) = 198.74 
but within la we still have useful information since these 
models are more probable than LCDM. 

We start with a simple generalization of LCDM: We 
replace the cosmological constant energy density by a 
dark energy with constant equation of state parameter. 
This ansatz has been called 'quiessence' in the literature 
[23. The form oi H{z) is 

H^Z; nora, W) = i/2[(f]„„(l + z)3 + (l_l]o„)(l+z)3(l+»)] 

(2.6) 

This ansatz has two free parameters flom and w. We 
use prior information from large scale structure (f2om^ = 
0.2 ± 0.03 m with h ^ 0.72 ± 0.08 IH) to fix = 0.3 
in this and in all subsequent ansatze. We thus evaluate 
X^(w) and minimize with respect to w. We find 

xLn = X'(^ = -0.93) - 199.3 (2.7) 

Including the la errors we have 

w = -0.93 ±0.08 (2.8) 

For f7o„i = 0.34 we findl23| 

xLn=x'(^« = -1-01) = 198.69 (2.9) 
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FIG. 2: The redshift dependence of the equation of state pa- 
rameter w{z) for the q — A ansatze. The thick curve is the 
best fit and the Ught (dark) shaded regions represent the la 
{2a) error regions. 



which is identical with the corresponding result of 
Ref . [isll . Thus, the minimization of this generalized 
ansatz gives a best fit that is indistinguishable at the 
Icr level from LCDM. This means either that LCDM is 
truly the best fit model or that we have not chosen a 
general enough ansatz to see a better fit. 

A further generalized ansatz involves the combination 
of cosmological constant with quiessence (quiessence-A 
ansatz). 

The form of H{z) in this case is 

H^iz; ai,wi) = H^[no-m{l + z)^ + ai(l + zf(^+-'^) + 
(1 - no,n - ai)] (2.10) 

Setting flom — 0.3 and minimizing x'^iO'i^'Wi) with re- 
spect to wi , ai we find 



= X^iwi = 3.44, ai = (5 lO^^) = 195.6 (2.11) 



Including the error bars we have 



ai ^ (5 l°-i) 10 



-4 



(2.12) 



Clearly the fit is better compared to LCDM but the 
X^{LCDM) = x^{ai = 0) = 200.3 corresponding to 
LCDM with r^om = 0.3 differs by less than Axl„{n 
2) = 6.17 from Xmin- Therefore, LCDM is consistent at 
the 2(7 level (but not at the la) with the best fit of this 
ansatz. Nevertheless, given that this fit is better it is 
interesting to compare the dark energy properties corre- 
sponding to this ansatz at best fit with those of LCDM. 
These properties are well described by the effective equa- 
tion of state parameter w(z) = £^ which in general (and 



in this case) depends on the redshift z. We can express 
w{z) in terms of H{z), ^ and fJom using the Friedman 
equations 



and 



SttG 



(Prn + PDe) 



(2.13) 



3i72 



+ {pDE + ipDE)] (2.14) 



where q is the deceleration parameter and we have de- 
fined as dark energy any other homogeneous and isotropic 
source of gravity apart from matter. Using (|2.13|) and 
(ITTl) we find 



PDE 



AtiG 

Using lfm| and (|TT^ we find 



, . ^ Pde{z) ^ 2q{z) - 1 
Pde{z) " 3(1 



where 



SnGpmjz) 
3iJ2(z) 



= rJom(l + z) 



3 -^0 



(2.15) 



(2.16) 



(2.17) 



Using now the definitions of q and H it is easy to show 
that 



-1 + il + z) 



dlnH 

dz 



Thus substituting (|2.18(l in H2.16|l we have 
Pde{z) _ |(l + z)-'^i"« ^ 



w{z) = 



dz 



PDe{z) 1 - (H^)2rtom{l + zY^ 



(2.18) 



(2.19) 



In the case of generalized Friedman equations valid in 
modified gravity models, equation l|2.19(l can still be use- 
ful in characterizing the expansion history but it should 
not be interpreted as a property of an energy substance. 
Using the best fit form of the quiessence-A (q-A) ansatz 
in (|2.19|l we find the predicted form of w{z) which is 
plotted in Fig. 2 along with the la and 2a error regions 
obtained by maximal variation of the parameters ai and 
wi within the la and 2a error contours of x^ as described 
in the previous section. This form of w{z) (without error 
regions) along with the corresponding forms predicted 
by the other ansatze discussed below, is also shown in 
Fig. 3. Clearly w{z) differs significantly from the LCDM 
prediction of w = — 1 at redshifts z > 0.5. In particu- 
lar we find w{z) ~ — 1 for z < 0.5 while ■w(z) ~ 3 for 

z ^ 1. Thus, this ansatz gives us a hint for the 'meta- 
morphosis' of dark energy from antigravity {w = —1) 
at low redshifts to 'hypergravity' {w ~ 3) at high red- 
shifts. Clearly this 'metamorphosis' (if true) can not 
persist to arbitrarily high redshifts due to constraints 
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coming from large scale structure and nucleosynthesis. 
Thus, it is either not realized in nature and we have 
w{z) < at all redshifts or it is part of an oscillating 
behavior of the dark energy equation of state param- 
eter. This later possibility could also help resolve the 
coincidence problem and is a prediction^ (see 
also '33,'3i|) of many models [H HE |H3 with stabilized 
modulus 38, 39] of extra dimensions. 

In addition to w{z) we also plot the reduced form of 
H{z) compared to LCDM defined as 



H\z) 



(2.20) 



where Hl^,j^j^j = 0.3 (1 + z)^ + 0.7. The reduced best 
fit H^{z) for the q-A ansatz is shown in Fig. 4 along 
with the best fits Hr{z) functions corresponding to some 
of the other ansatze discussed below. Even though both 
ansatze (q-A and LCDM) are consistent with the data 
and with each other at the 2a level the predicted forms 
of Hr{z) and 'w{z) differ significantly at z > 0.5. 

In order to shed more light to the dark energy 'meta- 
morphosis' puzzle we now consider more general forms of 
H{z) ansatze. For each ansatz we identify the predicted 
form of H{z) and the parameters requiring fitting. Then 
we evaluate and minimize with respect to these pa- 
rameters setting fiom = 0.3. Finally, we identify the best 
fit parameter values and the corresponding xinin (^^^ Ta- 
ble 1), plot the corresponding w{z) and H{z) (see Figs. 3 
and 4) and classify the models according to the goodness 
of fit. This classification will lead to some interesting 
conclusions about the generic properties of dark energy. 
The la and 2a regions of the w{z) and H{z) curves in the 
context of a particular ansatz are not particularly useful 
since other ansatze with equally good fits can give w{z) 
and H{z) best fits that are well outside the 2a regions of 
the initial ansatz especially in regions with z > 1 where 
just a few data points are available. Thus, in order to 
simplify the plots of Figs 3 and 4 and avoid confusion we 
only show the best fit curves without the corresponding 
la and 2a regions. 

We now briefly describe each one of the ansatze con- 
sidered: 

1. Cubic Polynomial in (1+z) (PS): 

H\z) = H^iQomil + zf + a3(l + zf + 02(1 + zf 
+ai{l + z) + {l-ai-a2- 03 - ^om)] (2.21) 

where ai, 02, 03 are unknown parameters to be 
fit. The only priors used in this (and the other 
ansatze discussed here) are fiom = 0.3 and flat- 
ness. In contrast to Ref. ^3 have not fixed 
aa = since large scale structure data do not ex- 
clude a non-clustering form of matter with prop- 
erties similar to those of hot dark matter with 
very large free streaming (or Jeans) length. This 
model has a slightly worse fit (Xmm = 196.6) com- 
pared to q-A (Xmin — 195.6) but it also corre- 



at z ~ 1. Its properties at best fit are shown in 
Figs. 3 and 4. The values of the best fit param- 
eters are oi = -2.55 ± 0.12, = 0.50 torn and 
as = 0.36±0.06. Compared to the quadratic ansatz 
of Ref. ^3 this cubic ansatz has slightly better 
Xmin and w{z) is larger than the corresponding 
w{z) of the quadratic ansatz at z > 1 shown in 
Ref. 40] and in Fig. 3. 

2. LinderflU] ansatz w{z) ^ wq + (LA): 

H^iz)^Hl[nomil + z)^+ 
(1 - nom){l + z)3(i+""'+"'i) e3™i(TT7-i)] (2.22) 

It can easily be verified using H2.19|l that the ansatz 
H2.22|l leads to a w(z) of the form 



w{z) = Wo + 



WiZ 

1 + z 



(2.23) 



This ansatz was suggested by Linder|4lJ and it is 
designed to interpolate between two values of w. 
w{z) = Wo at z ~ and w{z) ~ lui at z >> 1. It 
does not give a very good fit (compared to other 
two parameter ansatze) to the data {Xmin = 197.3) 
even though it is still better than LCDM. The rea- 
son is that it only allows a slow change of w{z) 
with redshift between wq and Wi while the data 
seem to require a more 'abrupt' change from wq 
to wi. The best fit values of these parameters are 
Wo = -1.29 ± 0.10 and wi = 2.84 ± 0.05. 

3. Linear Ansatz w{z) = wq + wiz: 

H\z) = Hl[no,n{l + zf + 

(1 - 17o™)(l + z)3(i+«'«-«'i) e^™^^] (2.24) 



This ansatz was suggested in Refs |42, |42l |4J] in a 
more general power series form. It can give very 
good fits for low redshift data (z 0.5) but it can 
not fit well the general form of w{z) at z ~ 1. At 
best fit it gives Xmm — 196.6 which is about av- 
erage compared to the other ansatze. The best 
fit parameter values are wq = —1.25 ±0.09 and 
wi = 1.97 torn which leads to w{z) > 0.5 for z ^ 1. 

4. Chaplyqinqas ( CG) and Generalized Ghaplyqiri gas 

/^CG^B El El 01 E IIS IM IM [M 111; Im IM 



H^z)=Hi[norn{l+z)^ + 



(1 - f^Om)V^+(l- A)(l+Z)^ 



(2.25) 



where the above form of H(z) is a generalization of 
the usual Chaplygin gas ansatz which is obtained 
for a = 6 0]. For a = 6 the equation of state of 
Chaplygin gas dark energy is 



sponds to w(z) 



-1 at z 0.5 and w{z) > 



Pc 



A 

Pc 



(2.26) 
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where 



FIG. 3: The redshift dependence of the equation of state pa- 
rameter for the cosmological ansatze of Table 1. The numbers 
in the parentheses indicate the value of Xmin for each ansatz 
and its rank according to increasing value of Xmm • The prior 
Oom = 0.3 has been used in all cases except LCDM where the 
best fit value Qom = 0.34 was used giving Xmin = 198.745. 



From the form of equation H2.25|l it is clear that 
the Chaplygin gas behaves like pressureless dust at 
high redshifts and like a cosmological constant at 
2: ~ 0. The sound velocity of Chaplygin gas grows 
rapidly and approaches the velocity of light at late 

redshifts {vs = = ^ - t^). Thus inho- 

mogeneities of Chaplygin gas do not grow (since 
the Jeans length approaches the horizon) and no 
constraints can be imposed from large scale struc- 
ture observations. The physical motivation of the 
Chaplygin gas equation of state (|2.2t)|) comes from 
string theories |57j|. In fact, considering a d brane 
in a 2 dimensional space-time, the introduction 
of light cone variables in the resulting Nambu-Goto 
action leads to the action of a Newtonian fluid with 
equation of state H2.26I) . 

The limitation of the Chaplygin gas ansatz is that 
it constrains w{z) to w{z) < at all flnite redshifts. 
Thus its goodness of flt is below average {Xmin — 
197.9 for A = 0.96 ± 0.03). This limitation does 
not exist for the generalized Chaplygin gas ansatz 
(arbitrary a) which gives a much better fit (x: 
195.2 for A = 0.9998 taoooi " = 17.68 
fit gives w{z) ~ 1.9 > for z ~ 1 (see Fig. 3) as do 
all the ansatze with above average goodness of fit. 

5. Generalized Cardassian Ansatz (CA): 

H^z) = H^iQomil + zf + (1 - nom).fxiz)] (2.27) 



■ - 
t'ofd- This 



fx{z) 



1-n 



^ 'Om 



1 



a/9 



Om (l + z)3(l-")9^ 

(2.28) 

This model |^ |E3 emerges from a generalization 
of the Friedman equations and predicts accelerated 
expansion at recent times without any dark energy. 
In this model the universe is flat and consists only 
of matter and radiation. Here we follow Ref. |^ 
and consider the generalization 12.281) of the orig- 
inal Cardassian ansatz of Ref. 10] . The original 
ansatz is obtained from equation (|2.28|l by setting 
q = 1 and is equivalent to quiessence {w = const.). 
The generalized Cardassian ansatz has been fitted 
to Snia data in Ref. 58] using a much smaller 
Snia dataset with redshifts z < 1. We find (see 
Table 1) that this ansatz gives a relatively poor fit 
to the data (x™„ = 198.6 for q = 0.025 1:q°o%, 
n — —23 tg) below the average goodness of fit for 
the models considered and only slightly better com- 
pared to LCDM with Qom = 0.3 (xLn = 200.4). 
The predicted w(z) increases with z but remains 
negative (see Fig. 3). 

6. Generalized Quardatic Ansatz (GQuad): 

H\z) = Hl[VLara{l + zf + ai(l + z)3(i+-0 + 
02(1 + z)3(i+-^) + (1 - ai - 02 - f^om)] (2.29) 

This is another generalization of the quadratic 
polynomial fit for H{z) of Ref. [l^]. Here we do 
not add an arbitrary cubic term. Instead we al- 
low the exponents of the two monomials to vary 
and minimize with respect to the four parameters 
fli, 02, wi and W2 instead of two parameters ai 
and 02 with fixed wi = —2/3 and W2 = ^1/3 
for the quadratic model. The fit of the generalized 
ansatz is better (Xmm = 195.2 for ai — 0.57 



02 = 0.003 



-1-0.0003 



0.03 
0.02' 



-0 00021 ^1 = —1.13 ± 0.24 and W2 — 
2.49 ± 0.02) than the quadratic ansatz (Xmin — 
196.8 for oi = -4.05 ti.lj and 02 = 1.79 toU) and 
the best fit form of w{z) differs significantly from 
the corresponding quadratic best fit particularly for 
z ~ 1. In particular we find that w{z) ~ —1 for 
z ~ 0.4 and w{z) ^ 2 for z ^ 1. For comparison 
the quadratic ansatz of Ref. j42|] predicts w{z) ~ 
for z ~ 1 with very small 2a errors at z > 1. 
This disagreement of our generalized ansatz with 
the quadratic ansatz at z > 1 despite the small 
2(7 error regions is another indication of the lim- 
ited usefulness of ploting la and 2a error regions 
of w{z) in the context of a particular ansatz. These 
regions can be easily violated in the context of an- 
other ansatz with better or similar fit. 

A variant of this ansatz is one where one of the 
two arbitrary power law terms is replaced by an 
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exponentially increasing term. The corresponding 
ansatz (Sinh+DE) is 

H^Z) = i/2[f^0rn(l + Zf + ai(l + 2)3(l+-i) + 

a2 sinh(w22:) + (1 - ai - ^om)] (2.30) 

The results for this ansatz are almost identical to 
those of equation (|2.29|) (see Table 1 and Fig. 3). 

7. Oscillating Ansatz (OA): 

H^{z) = H§[no„,{l + z f + ai cos{a2Z^ + 03) + 

(1 - aiCos(a3) - f^om)] (2.31) 

This is our best fit ansatz. It gives a better fit to 
the data than any of the other ansatze (Xmin = 

I 



Table 1 



Model 


H{z), (Oo™ = 0.3) 


■ 


Best Fit Parameters 


OA 
(1) 


H^{z) = H^[norn{l + zf + ai cos{a2z'^ + 03) 
+ (1 - ai cos(a3) - flom)] 


194.1 


03 = -0.067r ± O.OItt 


GCG (2) 


H-\z) = i/^[17o™(l + zf + (1 - no^)^A + (1 - ^)(1 + z)-] 


195.2 


A = 0.9998 ™,« = 17.68 t'of. 


GQ 

(3) 


H^z) = Hi[no,n{l + zf + ai(l + 2)3(i+-i) + 02(1 + z)3(i+-2) 

+ (1 - fli - a2 - flom)] 


195.2 


wi = -1.13 ± 0.24, W2 = 2.49 ± 0.02 
„ — nK7 +0 03 „ _ n nn-i+o oooa 
ai - u.57 _o.o2; a2 - 0.003_o 0002 


Sinh+DE 
(4) 


H'\z) = H'^[nn„Al + z)3 + ai(l + z)3(i+"'i) + a2 sinh(w2z) 
+1 - riom - ai] 


195.4 


= -0.81 +^-^, W2 = 5.90 +^•5'^ 
ai = -0.50 a2 = 0.026 t'oTil 


q-A (5) 


H^{z) - i/o'[f^o™(l + zf + ai(l + z)3(i+-i) + (1 - nora - ai)] 


195.6 


w, - 3 44 +°-'"' «i ~ 10-4 +0.4 10-* 


P3 
(6) 


H-'{z) = i/o'[f^Om(l + zf + 03(1 + zf + 02(1 + zf+ 
ai(l + z) + (1 - ai - 02 - 03 - r^om)] 


196.6 


ai = -2.55 ±0.12 
02 = 0.50 a3 = 0.36 ± 0.06 


Linear 
(7) 


H-'{z) = i/^ir!om(l + zf + (1 - 17o™) 

(^1 2;^3(l+wo-ii'i)g3mizJ 


196.6 


Wo = -1.25 ± 0.09, wi = 1.97 t^m 


Quad 
(8) 


H-\z) = i/^[r!o™(l + zf + ai(l + z) + 02(1 + zf 

+ (1 - fli - a2 - r^Om)] 


196.8 


a, = -4.05 tl'i?, a2 = 1-79 IJ]:^^ 


LA 
(9) 


H'\z) = i/o'[f^Om(l + zf + (1 - f^Om)(l + Z)3(l+-o+«.0 
g3ioi(^-l)j 


197.3 


Wo = -1-29 ± 0.10, wi = 2.84 ± 0.05 


CG (10) 


H-'{z) = i/o'[f^Om(l + ^)' + (1 - ^O.n)^A + (1 - A)(l + zf] 


197.9 


A = 0.96 ±0.03 


CA (11) 


H^{z) = i/o'[f^Om(l + zf + (1 - Oo„)/x(z)] 


198.6 


q = Om t^-Z, n = -23 tt 


LCDM (12) 


i/^(z) = H-^[nomil + zf + (1 - r!o™)J 


198.7 


rioni = 0.34 ± 0.03 


QUIES (13) 


H^Z) = Hi[norn{l + zf + (1 - Qomfl + z)3(l+-)] 


199.3 


w= -0.93 ±0.08 


SCDM 


iJ^(z) =iJ^(l+z)3 


431.4 





194.1 for ai = -3.36 Zq^^, 02 = 2.12 and 
03 = — 0.067r ± O.OItt). The behavior of w{z) how- 
ever for z 1 is qualitatively different compared 
to the other ansatze (see Fig. 3). For z 0.3 
we find w{z) ~ —1. For 0.5 ~ z ~ 1.2 we find 
wlz) > with a maximum w{z ~ 0.75) ~ 1.5. At 
z ^ 1.2, w{z) becomes negative and continues os- 
cillating around w ~ with large amplitude. This 
redshift range however includes only one data point 
at z = 1.75 which can not constrain the behavior of 
w(z) and H{z) in any statistically significant way. 



The quality of the fit for the oscillating ansatz com- 
bined with indications from the better fits of other 
ansatze that indicate w{z) > for z ~ 1 supports the 
idea that some type of oscillation probably takes place 
for w{z). 

From the theoretical viewpoint this idea is also sup- 
ported for two reasons: 

• Coincidence Problem Resolution: An oscillating ex- 
pansion rate can help resolve the coincidence prob- 
lem since our present accelerating phase is viewed 
simply as part of a sequence of accelerating and de- 
celeratin g p eriods in the expansion history of the 
universe |3lLl3^ 



• Extra Dimensions: Models with extra dimensions 
generically predict oscillations of the stabilized 
modulus of the extra dimension size (the radion 
field ) d ue to its coupling to redshifting matter 
prl H^ . These oscillations backreact on the ex- 
pansion rate and induce oscillations of the Hubble 
parameter. 

Another factor pointing towards oscillating expansion 
rate is the north-south pencil beam survey of Ref.|60j| 
which suggests an apparent periodicity in the galaxy 
distribution. The number of galaxies as a function of 
redshift seems to clump at regularly spaced intervals of 
128h~^Mpc. Recent simulations !61j have indicated that 
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0.25 0.5 0.75 1 1.25 1.5 1.75 

z 



FIG. 4: The reduced Hubble parameter for some of the best 
and the worst fits of the cosmological ansatze of Table 1. The 
number in the parenthesis shows the rank (1 — 13) of the 
corresponding ansatz in terms of goodness of fit. The LCDM 
curve is not flat at zero because in its construction we used 
the best fit value = 0.34 while the Hlcdm on the axis 
assumes the prior fiom = 0.30. 

this regularity has a priori probabihty less than 10^"^ in 
CDM universes with or without a cosmological constant. 
An oscillating expansion rate could resolve this puzzle 
without invoking special features in the primordial fluc- 
tuations spectrum. 



III. CONCLUSION 

We have fitted several cosmological models using the 
maximum likelihood method and the most recent Snia 
data consisting of 194 data points. No priors have been 
used in our analysis other than those indicated by other 
observations i.e. flatness and fiom = 0.3. The fact that 
we have fixed r^om instead of marginalizing over it could 
have artificially decreased somewhat the error bars of the 
parameters. This decrease however is not important in 
view of the main goal of this work which is the classifica- 
tion of the models studied according to their quality of 
fit to the data. 

We have confirmed recent studies indicat- 
ing an increase of the equation of state parameter with 
redshift termed 'metamorphosis' of dark energy in a re- 
cent study 40] . We have shown however that the best 
fit ansatze indicate that this 'metamorphosis' continues 
beyond w{z) — and leads to w{z) > at z ~ 1. Nucle- 
osynthesis and large scale structure constraints are not 
consistent with w{z) > at arbitrarily high redshifts. 
Thus our best fits to the data can only be made con- 



sistent with these constraints if some kind of oscillating 
behavior is realized for the effective dark energy equa- 
tion of state parameter w{z). This possibility is further 
enhanced by the fact that an oscillating expansion rate 
ansatz has provided the best fit to the data among all 
the 13 ansatze considered and also by other theoretical 
and observational arguments discussed in the previous 
section. 

At low redshifts (z < 0.5) all the fitted ansatze ap- 
proach w{z) ~ —1 with w{z) approximately constant. 
Most (but not all) of the better fits predict w{z) slightly 
less than —1 (up to —1.3) for some redshift range within 
[0, 0.5] but not necessarily at z = 0. There are good fits 
however (like the quiessence—A ansatz) with Xmin below 
average for which w{z) > —1 for all z implying that phan- 
tom energy ,6^] {w < —1) is consistent with the data 
but is not necessarily more probable than dark energy 
{w > —1). Since there are good fits with w{z) > —1 and 
rapidly increasing w(z) for z > 0.5, we conclude that even 
if the prior w > —1 were used, with the proper ansatz we 
would still be able to see the rapid increase of w(z) with z. 
This does not agree with the conclusion of Ref. [Iq] that 
it is the use of priors that would hide the increase of w{z) 
with redshift. Here we have shown that the cosmological 
ansatz selection can also play a crucial role in revealing or 
hiding the true expansion history of the universe. Other 
comparisons of particular ansatze with the SnIa data can 
also be found in the literaturefel IgtUsI IgIItoI . 

Thus we arc led to an important question: Is there 
a .systematic way to use the SnIa data in constructing a 
relatively simple cosmological ansatz that will give the best 
possible fit to the data given the number of parameters? 
Addressing this important issue will be the subject of a 
subsequent paper. 

Acknowledgements: We thank U. Alam , V. Sahni 
and J. Tonry for useful clarifications on the analysis of 
the SnIa data. This work was supported by the European 
Research and Training Network HPRN-CT-2000-00152. 

IV. APPENDIX 

Here we demonstrate that the marginalization over the 
zero point magnitude M defined in equation H1.6|l would 
have negligible effect (0(1%)) on our results. Any model 
will predict the theoretical value D]^{z;ai, ...,an) with 
some undetermined parameters {e.g. flrm^A)- The 
best- fit model is obtained by minimizing the auantitv|l8|: 

h (^^o,,oZ..(..))^ + ( "'°-"f;'"V ..)^ 

(4-1) 

where M' = M — Mobs is a free parameter representing 
the difference between the actual M, (see equation (|1.6() ). 
and it's assumed value Mobs in the data. 

Uniform marginalization over M' corresponds to inde- 
grating over M' and therefore working with a defined 
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by: 



with 



/+00 
-CO 

which after some manipulation gives: 

/ + 00 
-oo 

where 

N 2 ^ 



and 



^ 1 

C = 0.04^- 



(4.2) 



(4.3) 



(4.4) 



(4.5) 



Thus, the "marginahzed" over M' is: 



(4.6) 



x' = x'(M' = o) 



^ + lni^)^X^iM'^0) (4.7) 



because in the cases considered the last two terms on 
the RHS of equation (|4.7|) are of 0(1) while the first is 
O(IO^) and therefore dominates over the others. Thus 
the effects of the marginalization are of order 1% and 
can be neglected. This same conclusion has also been 
reached in Ref. 25]. 
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